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PREFACE 


This  Memorandum  is  a  result  of  RAND's  continuing  effort  to  develop 
improved  management  techniques  for  recoverable  items  (see  References). 
Past  RAND  work  on  stockage  policy  commonly  used  three  measures  of  sup¬ 
port  effectiveness  2  fill  rate,  backorders,  and  operational  rate.  Al¬ 
though  these  measures  are  fairly  well  understood  by  Air  Force  supply 
people  and  are  tractable  from  the  standpoint  of  computation,  it  is  dif¬ 
ficult  to  relate  them  directly  to  certain  measures  of  effectiveness  that 
are  more  operationally  meaningful. 

This  Memorandum  proposes  an  alternative  measure  of  effectiveness 
that  appears  to  be  a  better  predictor  of  the  expected  number  of  NORS 
aircraft  (Not  Operationally  Ready- -Supply) ,  and  provides  a  mathematical 
theorem  that  is  useful  in  minimizing  this  measure. 

Our  experience  in  using  this  new  measure  for  optimization  and  pre¬ 
diction  of  system  effectiveness  is  still  limited;  nevertheless,  initial 
results  indicate  that  the  characteristics  of  stockage  policies  obtained 
by  optimizing  the  NORS  criterion  have  a  great  deal  of  intuitive  appeal. 
Further  research  will  be  concentrated  on  developing  an  improved  algor¬ 
ithm  for  optimization. 

R.B.S.  Brooks  is  a  consultant  with  The  RAND  Corporation,  and  is 
affiliated  with  Bowdoin  College,  Maine. 


SUMMARY 


This  Memorandum  examines  four  measures  of  the  performance  of  the 
supply  system  at  an  air  base  in  providing  spare  parts  for  aircraft  re¬ 
pair:  fill  rate,  average  backorders,  operational  rate,  and  average  num¬ 
ber  of  NORS  aircraft. 

These  measures  are  briefly  defined  as  follows:  fill  rate  is  the 
ratio  of  the  number  of  units  Issued  over  a  fixed  time  period  to  the  num¬ 
ber  demanded  over  the  same  period.  Backorders  consist  of  the  number  of 
due-outs  from  base  supply  at  any  given  point  in  time.  Operational  rate 
Is  the  probability  that,  at  any  given  point  in  time,  there  will  be  no 
due-out  from  base  supply.  The  number  of  NORS  aircraft  is  the  number 
grounded  for  lack  of  spare  parts  at  any_ given  point  in  time. 

The  first  three  measures  have  often  been  used  in  RAND  studies  on 
stockage  policy  because  they  are  computationally  tractable;  and  the  Air 
Force  also  uses  fill  rate  as  a  yardstick  for  base  supply  performance. 

All  three  of  these  measures  suffer  a  common  drawback,  however:  it  is 
difficult  to  relate  them  to  operations. 

As  an  alternative  measure,  this  Memorandum  suggests  a  mathematical 
method  of  predicting  the  expected  number  of  NORS  aircraft.  Within  the 
limitations  imposed  by  assumptions,  this  NORS  criterion  brings  the  per¬ 
formance  myasure  for  base  supply  closer  to  the  measure  for  base  opera¬ 
tions.  Unfortunately,  the  NORS  criterion  is  computationally  cumbersome 
to  deal  with  because  it  Is  not  a  separable  function.  Hence,  item-by¬ 
item  optimization  with  a  cost  constraint  represented  by  a  Lagrange  multi¬ 
plier  is  not  feasible.  As  a  partial  solution  to  this  difficulty,  a 
theorem  that  shows  how  a  separable  function  may  be  substituted  for  expected 


NOH.S  and  an  algorithm  based  on  this  theorem  are  also  given.  Although 
experience  with  this  algorithm  is  still  limited,  it  seems  to  have  worked 
well  in  every  case  tried  at  RAND.  The  following  are  tentative  observa¬ 
tions  based  on  computational  experience! 

A.  The  NORS  criterion  function  is  almost  linear  for  a  wide 
range  of  cost.  This  means  that  a  small  change  in  the 
Lagrange  multiplier  representing  the  cost  of  a  NORS  air¬ 
craft  will  result  in  a  markedly  different  stoakage  policy. 

B .  In  terns  of  NORS,  the  NORS  policy  is  superior  to  the 
other  three  policies j  as  it  should  be.  The  operational 
rate  policy  seems  to  be  the  best  among  the  remaining 
three,  and  the  fill  rate  policy  is  the  least  satisfactory . 

C.  For  a  relatively  high  stockage  investment ,  the  four 
policies  yield  more  or  less  the  same  number  of  NORS . 

However,  when  the  constraint  on  investment  becomes  more 
stringent,  the  difference  in  the  characteristics  of  the 
four  policies  (measured  in  terms  of  the  predicted  NORS) 
becomes  more  pronounced,  i.e .,  the  performance  of  the 
three  non-NORS  policies  degrades  rather  rapidly  as  the 
investment  in  inventory  decreases.  This  degradation 
shows  up  before  investment  has  dropped  below  the  level 
needed  to  support  a  90  percent  fill  rate. 

D.  The  range  of  the  NORS  policy  is  less  than  those  of  the 
other  three  policies  for  a  relatively  low  investment ; 
however ,  as  investment  in  inventory  is  increased,  the 
NORS  policy  expands  its  range  more  quickly  than  do  the 
othere . 
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I.  INTRODUCTION 


This  Memorandum  compares  four  methods  of  measuring  the  supply  sys¬ 
tem  performance  at  an  Air  Force  base  in  providing  spare  parts  for  air¬ 
craft  repair.  The  analysis  will  be  limited  to  the  so-called  "recover¬ 
able"  items,  items  which  are  sufficiently  expensive  to  warrant  a  one-at- 
a-time  base  requisitioning  policy.  The  four  measures  are  fill  rate, 
average  backorders,  operational  rate,  and  the  average  number  of  NORS 
aircraft  (Not  Operationally  Ready--Supply  aircraft) . 

These  measures  are  defined  and  briefly  discussed  in  Sec.  II.  Sec¬ 
tion  III  reviews  the  necessary  inventory  theory  to  show  how  to  predict 
the  performance  of  the  supply  system  as  a  function  of  the  inventory  pol¬ 
icy  being  followed. 

Section  IV  examines  the  optimization  problem,  i.e.,  the  problem  of 
obtaining  a  policy  with  a  high  level  of  performance  at  a  low  cost.  This 
problem  has  a  fairly  well-known  and  simple  solution  in  the  case  of  the 
first  three  measures  of  performance.  Thus  far,  however,  it  has  been 
practically  impossible  to  optimize  the  NORS  measure  because  of  a  tech¬ 
nical  dif f lculty--lack  of  separability.  The  last  part  of  Sec.  IV  dis¬ 
cusses  this  difficulty  and  presents  a  partial  solution  to  it. 

Section  V  presents  a  numerical  comparison  of  policies  obtained  by 
optimizing  each  of  the  four  criteria.  Section  VI  gives  some  tentative 
conclusions.  The  proof  of  the  main  theorem  of  Sec.  IV  is  given  in  Ap¬ 
pendix  A.  Appendix  B  contains  a  description  and  listing  of  the  computer 
program  used  to  obtain  the  numerical  results  of  Sec.  V. 


XI.  MEASURES  OF  SUPPLY  PERFORMANCE 


This  section  describes  four  measures  of  supply  performance:  fill 
rate,  average  backorders,  operational  rate,  and  average  number  of  NORS 
aircraft. 

FILL  RATE 

Fill  rate  may  be  measured  by  taking  the  total  number  of  units  de¬ 
manded  from  base  supply  over  a  fixed  period  of  time --a  year,  say- -and 
dividing  that  number  into  the  total  number  of  units  issued  at  the  time 
they  were  demanded*  The  quotients  thus  the  percentage  of  demands  that 
were  immediately  filled.^ 

BACKORDERS 

A  backorder  is  defined  here  as  a  due-out  of  one  unit  of  stock  from 
base  supply  (not  a  due-in  from  depot).  To  measure  average  backorders, 
one  may  take  each  due -out  that  is  established  during  the  course  of  a 
year,  say,  observe  how  many  days  it  takes  to  satisfy  the  backorder,  add 
up  all  these  numbers  (one  number  for  each  due -out  established)  ,  and  di¬ 
vide  by  365.  Another  method  that  will  yield  almost  the  same  answer  is 
to  observe  the  number  of  due-outs  in  existence  at  a  fixed  time  each  day, 
and  then  average  all  fhese  numbers  together  over  the  course  of  e  year. 

Average  backorders  have  an  advantage  over  fill  rate  as  a  measure 
of  performance,  since  we  care  not  only  whether  backorders  occur,  but 
also  how  long  they  last*  To  take  an  extreme  example,  a  supply  system 
with  zero  fill  rate  will  stLll  be  very  good  if  each  backorder  lasts 

^In  this  RM  a  demand  is  a  request  for  issue  of  one  unit  of  stock* 
If,  in  a  single  request,  two  or  more  units  are  requested,  then  this  is 
counted  as  two  or  more  demands. 


only  three  minutes.  Fill  rate  gives,  In  this  case,  a  very  poor  indi¬ 
cation  of  performance.  On  the  other  hand,  since  the  average  number  of 
backorders  for  this  system  will  be  low  (unless  demand  rates  are  extra¬ 
ordinarily  high),  the  average  number  of  backorders  will,  in  this  case, 
be  a  good  measure  of  performance* 

OPERATIONAL  RATE 

Operational  rate  is  the  probability  that,  at  any  given  point  in 
time,  there  will  be  no  due-outs  from  base  supply  (backorders).  As  in 
the  case  of  average  backorders,  there  are  two  methods  for  measuring  op¬ 
erational  rate.  In  the  first  method  we  observe,  during  the  course  of  a 
year,  the  length  of  time  (in  days)  that  no  backorders  are  in  existence. 
That  is,  we  measure  (in  days)  the  length  of  each  interval  of  time  in 
which  there  are  no  backorders,  and  add  up  all  these  lengths  over  the 
course  of  a  year.  We  then  divide  this  number  by  365,  This  gives  us 
the  percentage  of  time  during  the  year  that  no  backorders  were  in  exis- 
tence--one  measure  of  operational  rate. 

The  second  method  is  to  observe,  at  a  fixed  time  each  day  during 
the  course  of  a  year,  whether  any  backorders  are  in  existence*  We  count 
up  the  number  of  days  on  which  no  backorders  were  observed,  and  divide 
by  365*  This  gives  us  another  measure  of  operational  rate,  which  should 
yield  approximately  the  same  results  as  those  under  the  first  method. 

Operational  rate  has  an  advantage  over  both  fill  rate  and  back¬ 
orders  in  that  it  may  be  directly  related  to  the  supply  system's  effect 
on  operations*  If  we  are  willing  to  assume  that  each  of  the  items  in 
the  set  of  Items  we  are  considering  is  essential  to  the  operation  of 


the  aircraft  ,  then  operational  rate  is  the  probability  that  no  aircraf  t 
will  be  lacking  an  essential  part-* in  other  words,  that  there  will  be 
no  NOR S  aircraft* 

Operational  rate  also  has  a  disadvantage  over  both  fill  rate  and 

* 

average  backorders  in  that  it  has  a  rather  bothersome  all-or-nothing 
character*  In  terms  of  the  NORS  interpretation,  operational  rate  does 
not  distinguish  among  the  situations  in  which  there  is  one  NORS  aircraft, 
two,  three,  and  so  on*  When  operational  rate  is  used  as  a  measure  of 
performance,  having  one  NORS  aircraft  is  as  bad  as  having* ten, 

A  variant--or  rather  a  whole  family  of  variants--of  operational 
rate  is  operational  rate  after  k  cannibalizations  .  where  k  is  some  whole 
number- -0 ,  1,  2,  etc.  Operational  rate  with  three  cannibal izat ions ,  say, 
is  measured  in  the  same  way  as  the  ordinary  operational  rate,  except 
that  now  we  allow  cannibal  izat  ions  of  up  to  three  aircraft  in  order  to 
satisfy  backorders.  This  has  the  same  effect  as  saying  that  the  parts 
on  up  to  three  aircraft  may  be  considered  as  belonging  to  the  spare  parts 
inventory . 

Just  like  ordinary  operational  rate,  operational  rate  after  k  can¬ 
nibal  izat  ions  has  an  Interpretation  In  terms  of  NORS*  We  may  think  of 
opera  tlonal  rate  after  k  cannibal  Izat  ions  as  the  probability  that  at 
any  given  point  in  time  there  will  be  no  more  than  k  NORS  aircraft, 

(This  is  under  the  assumption  that  parts  shortages  will  be  consolidated, 
by  means  of  cannibalization,  on  as  few  aircraft  as  possible*) 

AVERAGE  NORS  AIRCRAFT 

Average  NORS  aircraft,  just  like  average  backorders,  can  be  mea¬ 
sured  in  two  ways.  In  the  first  method  we  count,  for  each  aircraft, 


the  number  of  NGR5  days  during  the  course  of  a  year,  add  these  numbers 
up  (one  for  each  aircraft),  and  divide  by  365.  In  the  second  method 
we  observe  the  number  of  NORS  aircraft  each  day  at  a  specified  time. 

These  numbers  are  -then  averaged  over  the  course  of  a  year. 

Average  NORS  aircraft  is  almost  certainly  the  best  of  the  four 
measures  of  supply  performance  discussed  here,  for  after  all,  the  pur¬ 
pose  of  a  spare  parts  supply  system  is  to  maintain  the  operational  readi¬ 
ness  of  aircraft,  and  we  are  really  interested  in  the  performance  of  the 
system  only  insofar  as  it  does  or  does  not  achieve  this  purpose. 

The  main  disadvantages  of  average  NORS  aircraft  as  a  measure  of 
performance  have  to  do  with  mathematical  trac tability ,  A  mathematical 
prediction  of  average  NORS  aircraft  requires  more  restrictive  assump¬ 
tions  than  does  prediction  of  fill  rate,  average  backorders,  or  opera¬ 
tional  rate.  (These  assumptions  are  those  needed  to  fully  justify  the 
relation  pointed  out  above  between  operational  ready  rate  after  k  can- 
nlballzatlons  and  the  probability  of  no  more  than  k  NORS  aircraft*)  An 
even  more  bothersome  difficulty  is  the  purely  mathematical  problem  of 
optimizing  the  mathematically  predicted  average  number  of  NORS  aircraft. 
It  turns  out  that  a  purely  technical  difficulty  arises  here  that  may  be 
circumvented  when  optimizing  any  of  the  three  previous  measures.  Section 
IV  describes  this  difficulty  (lack  of  "separab ilityM)  and  presents  a  par¬ 


tial  solution. 


III.  MATHEMATICAL  PREDICTION  OF  SUPPLY  PERFORMANCE 


This  sections  describes  how,  on  the  basis  of  certain  assumptions, 
the  measures  of  supply  performance  discussed  in  Sec.  II  may  be  pre¬ 
dicted  for  any  given  base  stockage  policy.  Most  of  the  content  of 
this  section  is  a  review  of  the  inventory  theory  presented  in  earlier 
RAND  Metro randa  (see  Re ferences)  . 

The  section  first  discusses  the  prediction  of  individual  item  be¬ 
havior,  and  then  shows  how  these  predictions  may  be  used  to  predict  the 
measures  of  performance  defined  in  Sec.  II. 

PREDICTION  OF  INDIVIDUAL  ITEM  BEHAVIOR 

In  order  to  predict  the  behavior  of  an  individual  item,  we  will 
make  the  following  assumptions. 

Assumption  1:  One-for-One  Requisitioning 

Assume  that  whenever  one  unit  of  an  item  is  demanded  from  base 
supply,  either  a  replacement  for  it  is  requisitioned  from  the  depot, 
or  a  carcass  Is  turned  in  to  base  maintenance  for  eventual  repair  and 
return  in  a  serviceable  condition  to  base  supply.  In  short,  a  chain 
of  events  is  set  in  motion  that  will  result  in  the  eventual  replacement 
of  the  demanded  item.  If  the  demand  is  filled  from  stock  on  hand,  then 
the  replacement  will  be  used  to  replenish  stock  on  hand.  If  the  demand 
results  in  a  backorder  (due-out  from  base  supply)  ,  then  the  replacement 


will  be  used  to  satisfy  that  due-out. 
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Assumption  2;  Backordering  of  Unsatisfied  Demands 

If  a  demand  occurs,  and  a  unit  of  the  item  is  in  base  supply,  that 
unit  will  be  used  to  satisfy  the  demand.  Otherwise,  the  demand  results 
in  a  backorder  or  due-out  from  base  supply. 

Assumption  3:  Markov  Property  for  Demand 

The  number  of  demands  for  an  item  that  occur  within  any  given  in¬ 
terval  of  time  is  a  random  variable;  this  number  is  statistically  inde¬ 
pendent  of  the  number  that  occur  in  another  interval  of  time,  unless  the 
two  intervals  overlap. 

Assumption  4:  Stationarity  of  Demand 

The  number  of  demands  within  an  interval  of  time  is  a  random  vari¬ 
able  whose  probability  distribution  depends  only  on  the  length  of  the 
interval.  Thus,  the  number  of  demands  during  January  has  the  same  prob¬ 
ability  distribution  as  the  number  in  December. 

Assumption  5:  Independence  of  Resupply  Time  and  Demand 

We  assume  that  the  length  of  time  required  for  a  requisition  on 
the  depot  to  be  filled  is  a  random  variable  independent  of  the  number 
of  demands  that  occur  in  any  interval  of  time,  and  similarly  for  the 
length  of  time  required  to  repair  an  item  on  the  base.  Finally,  the 
decision  on  whether  to  repair  an  item  on  base  is  independent  of  the  num¬ 
ber  of  demands  that  occur  in  any  given  interval  of  time. 

Under  the  above  assumptions  about  demand  (Assumptions  3  and  4)  ,  it 
can  be  shown  [14]  that  the  number  of  demands  that  occur  in  any  given  in¬ 
terval  of  time  has  a  compound  Poisson  probability  distribution 


e  ^ ,  for  k  =>  0 , 


f(k;t)  = 


E  5 ^  g(*n\k)  ,  for  k  >  0 

10  . 

_ 


where  t  is  the  length  of  the  time  interval ,  X  is  a  constant,  g  is  a  dis- 

(*n> 

Crete  probability  distribution  function,  and  g  denotes  its  n-£oid  con* 
volution. 

Under  our  assumptions  about  requisitioning  and  backordering  (assump¬ 
tions  1  and  2) ,  the  number  of  units  on  hand  plus  the  number  that  are  due 
in  (either  from  depot  or  base  maintenance)  minus  the  number  of  backorders 
remains  the  same.  This  constant  is  called  the  level  of  the  item: 


on  hand  4-  [due- ins]  -  backorders  =  level 


Moreover,  either  the  quantity  on  hand  is  zero,  or  the  number  of 
backorders  is  zero.  Thus,  if  we  know  the  level  and  also  know  the  due- 
ins,  we  also  know  both  on-hand  assets  and  backorders.  Explicitly: 


On  hand  = 


levei-[due-lns] ,  if  level  >  due-ins 
0,  if  level  <  [due- ins] 


Backo rders 


[ due -ins]- level ,  if  [due- ins]  >  level 
0,  if  [due- ins]  <  level 


Hence,  we  completely  know  the  state  of  the  item  if  we  know  its  level, 
which  is  a  matter  of  policy,  and  its  due-ins,  which  is  a  matter  of 


chance.  The  remarkable  fact  la  that  at  any  point  in  time,  under  the 
assumptions  we  have  made,  the  quantity  of  due-ins  is  a  random  variable 
whose  probability  distribution  is  given  by 

P{ [due- ins]  »  k]  *  f(k;T) 


where  T  is  the  average  resupply  time  for  the  item.  That  is,  T  is  the 
average  requisitioning  time  for  an  item  that  is  never  repaired  on  base; 

T  is  the  average  base  repair  time  for  an  item  that  is  always  repaired 
on  base;  and  T  is  a  weighted  average  of  these  two  times  for  an  item  that 
is  only  sometimes  repaired  on  base — where  the  weight  accorded  average 
base  repair  time  is  the  percentage  base  repair.  This  fact  is  the  ex¬ 
tended  form  of  Palm's  theorem  given  in  [6]. 


PREDICTION  OF  MEASURES  OF  PERFORMANCE 

We  now  suppose  that  we  have  n  reparable  items  on  the  base,  each  of 


which  obeys  the  foregoing  assumptions.  We  let  f^  be  the  probability  dis 

tribution  for  the  number  of  demands  for  the  jth  item  during  a  length  of 

time  equal  to  the  average  resupply  time  of  the  item.  Thus,  f^(k)  ts> 
til 

for  the  j  item,  the  probability  f(k;T)  introduced  above.  We  will  also 
let  be  the  level  for  the  j  item  and  d^  be  the  average  annual  demand 


Fill  Rate 

The  on-hand  assets  of  item  j  at  time  t  will  depend  only  on  demands 
previous  to  time  t,  and  not  upon  demands  at  time  t  or  afterward.  Thus 
the  proportion  of  demands  that  occur  when  the  on-hand  assets  are  zero 
will  be  the  same  as  the  proportion  of  time  that  on-hand  assets  are  zero. 
Hence,  the  proportion  of  demands  that  cannot  be  met  from  on-hand  assets 


is  equal  to  the  proportion  of  time  that  on-hand  assets  are  zero,  Or, 
to  put  the  matter  the  other  way  around,  the  proportion  of  those  demands 
for  item  j  that  can  be  filled  from  on-hand  assets  Is  equal  to  the  pro¬ 
portion  of  time  that  on-hand  assets  are  positive.  But  this  proportion 
is  just  the  probability  that  at  a  given  point  in  time  on-hand  assets  are 
positive.  As  pointed  out  in  the  previous  subsection,  on-hand  assets  will 
be  positive  if  and  only  if  the  number  of  due- ins  is  less  than  the  level. 
By  the  generalized  form  of  Faints  theorem  referred  to  above,  this  prob¬ 
ability  is  simply 


E  £  (k). 


We  will  let  Fj  be  the  cumulative  probability  distribution  corresponding 
to  f^  j  so  the  above  discussion  may  be  summarized  by  saying  that  the  fill 
rate  for  item  j,  i.e.,  the  proportion  of  demands  for  item  j  that  are 
filled  from  on-hand  assets,  is 


r.Cq.-l)  =  S  £.(k). 
J  J  '-'j J 


This  being  the  case,  the  average  number  of  demands  for  item  j  that  will 
be  filled  from  on-hand  assets  during  the  course  of  a  year  will  be  the 
above  proportion  multiplied  by  the  average  number  of  demands  during  the 


year.  Thus, 


is  the  average  number  of  demands  for  item  j  that  will  be  filled  during 
the  course  of  a  year.  When  we  add  these  numbers  up  over  all  the  items 
we  obtain 

3iwvi> 

as  the  total  number  of  demands  for  all  Items  that  are  filled  from  on- 
hand  assets  during  the  course  of  a  year.  In  order  to  obtain  the  fill 
rate  over  all  items  *  we  simply  divide  this  number  by  the  average  number 
of  demands  across  all  items: 

a  /  n 

Fill  rate  =*  E  d  F.(q  -1)/  E  d  . 

j«l  J  J  J  /  j-1  3 


Average  Backorders 

From  the  previous  subsection,  the  probability  of  k  backorders, 
k  >  0,  for  item  j  is  the  same  as  the  probability  of  k  +  due-ins  for 
item  j,  and  this  probability  is  simply 

Thus,  the  average  number  of  backorders  for  item  j  is 

to 

E  f(k+q  )k. 
k=l  3 


In  order  to  obtain  average  backorders  for  all  items,  we  simply  add  up 
the  average  backorders  for  the  individual  items. 


Average  Backorders 


n  « 

£  £  f(k+q.)k. 

j-1  k=l  J 


Operational  Rate 

In  order  to  predict  operational  rate  we  shall  have  to  make  one  more 
assumption  in  addition  to  those  we  have  made  already: 

Assumption  6;  Independence  of  demand  and  resupply  times  among  items. 
We  assume  that  the  number  of  demands  for  each  of  the  items  during  any 
fixed  interval  of  time  forms  a  family  of  independent  random  variables. 

We  also  assume  that  the  resupply  times  from  one  demand  to  another  are 
independent  when  the  demands  are  for  different  items. 

Given  this  assumption,  the  number  of  due- ins  for  the  various  items 
form  a  family  of  independent  random  variables.  Hence  the  probability 
that  there  will  be  no  backorders  at  all  is  the  product  over  all  the  items 
j  of  the  probability  of  no  backorders  on  item  j.  The  probability  of  no 
backorders  on  item  j  is  the  probability  that  the  number  of  due -ins  for 
item  j  is  less  than  or  equal  to  Its  level.  This  probability  is  simply 
F^(qj).  Hence, 

n 

Operational  rate  «  II  F  (q  )  , 

J-i  3  3 

When  one  aircraft  is  available  for  cannibalization,  it  is,  for  com¬ 
putational  purposes,  as  though  an  aircraft  had  been  disassembled  and  its 
parts  added  to  the  spare  parts  stock  level.  Let  a^  be  the  number  of  ap¬ 
plications  of  item  j  on  a  single  aircrafts  Then,  when  one  aircraft  la 


made  available  for  cannibalization,  the  stock  level  for  item  j  has  been, 


in  effect,  Increased  from  q. 


to  +  .J 


Thus  , 


Operational  rate  given  one  n 
aircraft  available  for  =  n 
cannibalization  j*l  ^  ^ 


In  general,  if  there  are  k  aircraft  available  for  cannibalization, 
then  the  stock  level  for  item  j  is,  in  effect,  Increased  from  q^  to 
q^  +  ka ^ ,  Thus , 


Operational  rate  given  n 

k  aircraft  available  =  n  F,(q,+ka  ). 
for  cannibalization  j=l  J  ^ 

Average  NORS  Aircraft 

In  order  to  predict  average  NORS  aircraft  we  need  two  more  assump¬ 
tions  In  addition  to  those  needed  up  to  this  point.  They  are: 

Item  essentiality.  Each  of  the  items  in  the  range  of  items  ia  es¬ 
sential  to  the  operation  of  the  aircraft;  and  any  item  that  could  cause 
a  NORS  is  in  this  set. 

Consolidation  of  parts  shortages  (cannibalization) .  When  a  demand 
for  a  part  cannot  be  satisfied  from  stock  on  hand,  then  it  is  satisfied, 
if  possible,  by  removing  the  needed  part  from  an  already  NORS  aircraft. 

In  fact,  parts  shortages  are  consolidated  on  as  few  aircraft  as  possible. 

Under  this  last  assumption,  if  we  start  from  a  situation  In  which 
there  are  no  NORS,  but  then  a  parts  shortage  occurs,  the  resulting  NORS 
aircraft  may  be  used  to  satisfy  other  parts  shortages  until  the  number 
of  backorders  on  one  of  the  parts,  part  j  say,  exceeds  a  at  which  point 


a  second  aircraft  will  go  NORS.  Thus,  the  probability  of  one  NORS  air¬ 
craft  or  fewer  is  simply  the  probability  that  for  all  parts  j  the  num¬ 


ber  of  backorders  on  part  j  does  not  exceed  a ^ .  Since  the  number  of 

backorders  on  part  j  is  0  unless  due-ins  exceed  and  is  (due-ins) -q^ 

otherwise,  the  probability  that  there  will  be  no  more  than  backorders 

on  part  j  is  the  probability  that  the  number  of  due- ins  on  part  j  does 

not  exceed  q.  +  a. ,  and  this  probability  is  simply  F.(q.  +  a.).  Thus 
J  J  J  J  3 

Fj(qj  +  a^)  is  the  probability  that  there  will  be  no  more  than  a^  back¬ 
orders  on  part  j.  The  probability  that  this  will  be  true  for  all  parts 
j  is  obtained  by  taking  the  product  of  these  probabilities.  Thus, 


Probability  of  one  or  n  F  fa  +a  V 
less  NORS  aircraft  “  ,pHj Y 

33^ 


This  is  the  same  number  as  the  operational  rate  given  one  aircraft  is 
available  for  cannibalization.  The  point  is  that  there  is  no  need  to 
cannibalize  an  aircraft  unless  there  is  a  parts  shortage,  so  the  prob¬ 
ability  that  there  is  at  most  one  NORS  is  the  same  as  the  probability 
that  all  demands  can  be  met  either  from  stock  on  hand  or  by  creating 
holes  in  at  most  one  aircraft,  and  this  is  simply  the  operational  rate 
given  one  aircraft  is  available  for  cannibalization. 

In  the  same  way,  we  have  the  general  result  that  the  probability 
that  there  are  k  or  less  NORS  is  the  same  as  the  operational  rate  given 
k  aircraft  are  available  for  cannibalization.  Thus, 

n 

-  HF(q+ka>. 

J-l  J  J  J 


Probability  of  k  or 
less  NORS  aircraft 


This  last  equation  gives  us,  then,  the  cumulative  probability  dis- 


ribution  function  for  the  number  of  NORS  aircraft.  Now  it  is  a  general 
fact  that  if  G  is  the  cumulative  probability  distribution  function  for 
a  nonnegative  discrete  random  variable  X,  then  the  expectation  of  X  may 
be  obtained  by  summing  the  ‘’tail11  of  G,  i.e., 

GD 

E(X)  »  Z  (l-G(k)). 
k=0 

This  fact  together  with  the  previous  equation  gives  us 

<*>  n 

Expected  NORS  =  £  (1-  II  F.(q.+ka.)). 

k=0  j=»l  J  J  J 


IV  *  OPTIMIZATION  OF  SUPPLY  PERFORMANCE 


In  this  section  we  discuss  the  problem  of  finding  stockage  poll- 
cies  that  optimize  supply  performance  and  cost.  A  policy  is  described 

by  specifying  the  levels  . q^  of  each  of  the  n  items.  The  cost  of 

a  policy  (q^. .  .  ,qn)  is 

n 


The  performance  is  measured  by  one  of  the  four  measures  described  pre¬ 
viously:  fill  rate,  backorders,  operational  rate,  or  NORS, 

EFFICIENT  POLICIES  AND  LAGRANGE  MULTIPLIERS 

To  facilitate  the  following  discussion,  let  Q  be  the  set  of  all 
possible  stockage  policies,  so  Q  Is  the  set  of  all  n- tuples  q  *  (q  , 
q^>  of  nonnegative  integers  We  write  C(q)  for  the  cost  of 

a  particular  policy  q.  Thus, 

c(q! . v  -  !cjY 

We  write  F(q)  for  the  performance  of  a  policy  q.  We  will  assume  that 
high  values  of  P  are  good  whereas  low  values  are  bad,  (In  the  case  of 
backorders  and  NORS  we  may  replace  F  by  -P,)  A  policy  q  is  said  to  be 
efficient  if  there  is  no  other  policy  which  is  better  on  the  basis  of 
performance,  and  no  worse  on  the  basis  of  cost,  and  also  there  Is  no 
policy  which  Is  better  on  the  basis  of  cost  and  no  worse  on  the  basis 
of  performance.  More  formally,  q  is  efficient  if,  given  any  other 


policy  q',  the  inequalities  C(q)  >  C(q')  and  P(q)  <  P(q')  together  Imply 
C(q)  *  C(q')  and  P(q)  *  P(q'). 

Efficiency  is  certainly  something  to  be  desired  of  a  policy,  for 
if  a  policy  is  not  efficient,  then  it  can  be  improved  upon  according  to 
at  least  one  criterion  (cost  or  performance)  without  degradation  with 
respect  to  the  other.  The  theorem  below  goes  back  in  one  form  or  ano¬ 
ther  to  the  eighteenth-century  mathematician  Joseph  Louis  Lagrange,  and 
is  given  in  the  following  form  in  Everett's  article  [3]: 

THEOREM  1.  Suppose  X  >  0,  and  q*  €  Q  maximizes  P(q)  -  XC(q)  , 
q  €  Q.  Then  q*  ^s  efficient. 

The  proof  is  quite  simple:  Suppose  q*  were  not  efficient.  Then 
there  would  be  q*  €  Q  with  P(q')  >  P(q*)  ,  C(q')  <  C(q*)  ,  and  at  least 
one  of  the  inequalities  would  be  strict.  For  this  q'  we  would  then  have 
P(q')  -  XC(q')  >  P(q*)  -  XC(q*)  ,  contrary  to  hypothesis.  Thus  q*  must 
be  efficient. 

The  number  X  in  the  above  theorem  is  called  a  "Lagrange  Multiplier." 
We  are  not  solely  interested  in  finding  an  efficient  policy,  however. 

For  example,  the  policy  in  which  all  item  levels  are  zero  is  efficient, 
since  every  other  policy  costs  more,  but  in  most  instances  we  would  not 
be  satisfied  with  the  performance  of  such  a  policy.  Actually,  we  are 
interested  in  finding  efficient  policies  that  will  meet,  at  least  ap¬ 
proximately,  a  prespecified  cost  or  a  pre specified  performance.  It  is 
possible  to  find  such  policies  by  varying  the  value  of  X  in  the  above 
theorem;  high  values  of  X  will  result  in  low-cost  low- performance  poli¬ 
cies;  low  values  of  X  will  result  in  high-cost  and  high-performance  poli¬ 
cies.  Systematic  methods  of  varying  X  so  as  to  meet  a  pre spec if  led  cost 


or  performance  level  may  be  found  in  [4],  [93,  [ll  and  [10]*  The  first 
two  references  use  what  has  sometimes  been  called  incremental  allocation 
or  marginal  analysis.  The  third  reference  shows  how  linear  programming 
may  be  used  to  find  the  right  X.  The  last  reference  suggests  a  binary 
search  technique.  Here,  we  are  interested  in  the  problem  of  optimizing 
P(q)  ~  XC( q)  for  fixed  X. 

THE  ROLE  OF  SEPARABILITY 

A  function  ¥  of  several  variables  is  called  separable  if 

there  are  n  functions  ¥-*¥„,*,,.¥  of  one  variable,  such  that 

1 J  2  1  n 

. v  -  jiw- 

A  very  pleasant  property  of  separable  functions  is  that  they  may  be  maxi¬ 
mized  one  variable  at  a  time*  That  is,  if  ¥^  attains  its  maximum  at  X*, 
^2  attains  its  maximum  at  x*,  and  so  on,  then  Y  attains  its  maximum  at 
(x*  ,x*  ,  * , .  ,x^)  ,  This  property  can  lead  to  considerable  savings  in  com¬ 
putational  time;  in  fact  it  can  make  the  difference  between  a  possible 
and  an  impossible  computation.  For  example,  if  n  =  1000  and  each  of  the 
variables  x^ > * ■ * >xn  can  on  10  different  values,  then  we  need  only 

look  at  10  values  of  in  order  to  maximize  4^,  10  values  of  x^  in  order 
to  maximize  ,  find  so  on,  leading  all  in  all  to  10,000  different  values 
of  the  variables  .  sx  .  This  is  a  possible  task.  If  ¥  is  not  sep¬ 

arable,  however,  we  must  look  at  all  the  10^^  different  values  that  the 
n- tuple  (x^,.,.,x^)  can  take  on.  This  is  an  impossible  task. 


It  Is  clear  that  if  two  functions  are  separable,  then  so  is  their 
sum  and  difference.  In  fact,  any  linear  combination  of  separable  func¬ 
tions  is  separable.  It  is  also  clear  that  the  cost  C(q)  =  2c, q,  of  a 

j  2  3 

stockage  policy  q  =  (q^,...,qn)  is  a  separable  function  of  the  stock 
levels  qt,...,qn.  Hence  P(q)  -  XC(q)  will  be  a  separable  function  of  the 
stock  levels  if  P,  the  measure  of  performance,  is  separable. 

In  the  case  of  fill  rate  we  have: 

n  n 

PCV-.-.O  «  2  d  F  (q  -1)/  2  d 
1  j»l  J  J  J  j=*l  J 

which  is  clearly  a  separable  function  of  the  individual  item  levels 
•  ,qn .  In  the  case  of  backorders 

n  « 

j  ■  *  *  jO  —  ^  E  E  f  (k-Ki  }  k , 

1  j=*l  k=l  J 

which  is  also  separable-  Now  operational  rate 

n 

n  F  (q  ) 
j=l  J  J 

Is  not  a  separable  function,  but  this  problem  is  easily  overcome  by 
working  with  the  logarithm  of  operational  rate,  instead  of  operational 
rate  itself;  for  certainly,  the  bigger  operational  rate  is,  the  bigger 
its  logarithm  is,  and  conversely.  Hence,  to  optimize  the  logarithm  of 
operational  rate  is  the  same  as  to  optimize  operational  rate.  The  log¬ 


arithm  of  operational  rate  is 


More  generally,  the  logarithm  of  operational  rate  given  k  aircraft  avail¬ 
able  for  cannibalization  is 

n 

P(q  ,...,q  )  =  E  log  F(q ,+ka.)  , 

J-l  3  3 

and  these  are  clearly  separable  functions  of  the  stock  levels.  Thus, 
Lagrange  multiplier  techniques  coupled  with  separability  of  the  perfor¬ 
mance  criteria  enable  us  to  solve  the  optimization  problem  in  the  case 
of  fill  rate,  backorders,  and  operational  rate. 

In  the  case  of  expected  NORS  aircraft,  however,  we  have 

“  n 

V(q  =  -  E  Cl-  n  F  (q  +ka  )), 

k=0  j=l  J  3  J 

* 

which  is  not  a  separable  function  of  the  item  levels.  Nor  can  we  use 
the  same  trick  used  on  operational  rate  in  order  to  turn  NORS  into  a 
separable  function.  This  lack  of  separability  is  the  main  reason  that 
NORS  cannot  be  used  for  a  performance  criterion  in  selecting  an  optimal 
policy.  The  next  subsection  presents  a  partial,  but  by  no  means  complete, 
solution  to  the  problem  of  the  lack  of  separability. 

OPTIMIZING  NORS 

We  show  here  that  in  optimizing  NORS,  it  is  possible  to  replace  the 


NORS  function 


2  (1-1IF  (q.-Hsa  )) 
k-G  j  J  J  J 


by  a  linear  combination  of  separable  functions.  Unfortunately *  there 
seems  to  be  no  systematic  method  for  finding  the  coefficients  in  the 
linear  combination. 

We  first  note  that  since  the  above  sum  is  supposed  to  converge,  we 
may  replace  it  by  a  finite  sum 

K 

E  (1-flF  (q .+ka  ))  , 
k=0  j  J  J  J 


and  the  mistake  committed  by  so  doing  may  be  made  arbitrarily  small  by 
making  K  sufficiently  large.  But  this  last  sum  is  simply 


K 

K-  Z  nF(q+kaJ, 
k=0  j  J  J  J 


Since  the  constant  K  is  unaffected  by  the  stockage  policy,  we  may  use 


K 

Z  IIF  (q  +ka.) 
W  j  J  J  J 


as  a  measure  of  performance. 

Because  we  are  taking  a  product  rather  than  a  sum*  our  measure  of 
performance  is  nonseparable «  If*  however,  we  could  somehow  replace  the 
term  IlFjCq^+ka^}  by  some  multiple  of  its  logarithm*  then  we  would  have 
a  separable  function ♦  We  are  thus  led  to  Che  following  question: 


Are  there  numbers  b^,.. 

K 


,bR  such  that  any  that  maximizes 


£  log  F  (q  +a  k)  -  X£c  q 
j=l  J  J  J  j  J  J 


also  maximizes 


K  n  n 

£  n  F  (q.+a  k)  -  x  2  C  q  1 
k=0  j=l  J  3  J  j-1  •»  J 


The  answer  to  this  question  is  yes,  as  the  following  theorem  and  its 
corollary  establish. 

THEOREM  2 .  Suppose  q* , . . . ,q*  are  nonnegative  Integers,  and  define 

numbers  .  ,b„  by 

—  U  K 


b  =  n  F . (q^+ka ,) . 

*  j-1  J  J  3 


Then  for  any  other  nonneeative  integers  q^ , . . . ,qn  we  have 

K  n  n  K  n  n 

£  n  F.(q.+ka  )  -  X  £  c  q  >  £  II  F.(q*+ka  )  -  X  £  c  q* 

k=0  j»l  J  J  3  j=l  33  k-0  j-1  J  J  J  j«l  J  J 


whenever 


K  n 

£  b  £  log  F  (q  +ka  ) 
k=0  j-1  J  J  J 


n 


K  n 


x  2  c.q  >  £  b  £  log  F(q*+ka  )  -  X  £  c  q*  . 
j-1  J  J  k-0  Kj-1  3  3  j*l  J  J 


In  particular,  when  q^j..«,q*  actually  optimizes  NORS  we  have  the 


answer  to  our  question. 


COROLLARY ,  Suppose  that 


K  n  n 

E  n  F,(q,  +ka . )  -  X  2  c.q. 

k-0  j-1  J  ■>  J  j.l  J  J 


actually  attains  a  maximum  when  (q  , , . , ,q^)  -  ( q*  , , . *  ,q*) *  Then  there 
are  nannegative  numbers  b^  ,  -  ,  *  ,b^  such  that  any  q**  f , ,q£*  that  maxi¬ 
mizes 


JT  ^  tl 

E  b  2  log  F  (q  +ka  )  -  X  E  c  q 
k=0  j=l  21  2  j-1  J  3 


also  maximizes 


K  n  n 

£  II  F  (q  +ka  )  -  X  £  c  q  , 
lc=0  j^X  J  J  J  J  J 


Xn  fact  the  numbers  b  *  * ,  *  #b  may  be  defined  as  in  the  theorem. 

"  '  U  Jv 

Thus,  instead  of  optimizing  the  nonseparable  function 


K 

£  nF  (q  +a  k)  , 
k=0j  J  J  J 


we  may  achieve  the  same  result  by  optimizing  a  separable  function 


K  n 

£  b  E  log  F  (q  +a  k>  . 
k=0  Kj=l  J  J  J 


Unfortunately,  in  order  to  do  this  we  need  to  know  the  values  of  . 

b  .  Xn  order  to  apply  our  corollary  to  determine  b  we  need  to 

U  lx 


know  the  numbe  rs 


-Z0“ 


for  various  values  of  \,  until  the  prespecified  cost  is  again  met.  We 

then  reset  bQ . and  repeat  the  process.  These  steps  are  repeated 

until  the  values  of  bg  , . .  .  ,b  stabilize.  The  adjustment  of  \  was  done 
by  linear  prog  ramming- -as  suggested  in  [l]. 

For  each  case  we  have  tried,  we  have  started  with  two  extreme  val¬ 
ues  for  (bg,...,b  ),  The  first  extreme  value  is  given  by 


b 


0 


K-l 


0,  b. 


1. 


This  corresponds  to  the  very  pessimistic  view  that  most  demands  will  be 
satisfied  by  cannibalization  of  NORS  aircraft.  The  resulting  policy 
(until  the  numbers  bn,...,b  are  subsequently  modified)  has  a  limited 
range  of  items  with  positive  levels  but  greater  depth  on  high-demand 
items.  The  other  extreme  set  of  values  is 


This  corresponds  to  the  case  in  which  we  think  there  will  be  virtually 
no  NORS  aircraft,  so  all  demands  will  have  to  be  satisfied  from  stock 
rather  than  cannibalization.  The  resulting  policy  (until  the  numbers 
b_,.,.,b  are  modified)  has  a  large  range  of  items  with  positive  levels, 
at  the  expense  of  depth  on  high-demand  items. 

In  all  the  cases  we  have  tried,  we  have  converged  to  the  same  val¬ 
ues  of  bg,...,bR  when  we  started  with  the  pessimistic  extreme  as  we  did 
in  starting  with  the  optimistic  extreme. 


V ,  NUMERICAL  RESULTS 


The  previous  sections  discussed  four  measures  of  supply  performance 
and  indicated  that  the  NORS  criterion  is  the  roost  desirable  because  it 
is  more  operationally  oriented  than  the  others.  To  give  some  quantita- 
tive  indication  of  how  the  stockage  policies  based  on  these  four  mea¬ 
sures  differ  from  each  other,  this  section  presents  some  numerical  re- 
sul ts . 

For  varying  investment  levels,  we  computed  four  stockage  policies, 
each  of  which  optimizes  one  of  the  four  criteria.  We  then  evaluated  the 
fill  rate,  backorders,  and  operational  rate  policies  for  their  respec¬ 
tive  NORS.  We  also  computed,  for  each  policy,  the  proportion  of  items 
that  were  given  a  positive  stock  level.  We  call  this  characteristic  the 
range  of  a  stockage  policy.  These  results  are  given  in  Table  2.  Before 
we  interpret  the  results  a  brief  description  of  the  data  used  for  the 
computations  is  in  order. 

DATA 

Table  1  summarizes  the  item  data,  obtained  by  aggregating  the  F-101 
data  we  used  for  a  field  test  of  the  RAND  Base  Stockage  Model  [5]  at 
Hamilton  Air  Force  Base.  The  488  items  are  grouped  into  189  different 
"item  types."  The  items  within  each  type  all  have  approximately  the 
same  unit  cost,  demand  rate,  and  response  time  (resupply  time).  The 
second  column  shows  the  number  of  items  in  each  item  type.  The  third 
column  shows  the  average  unit  cost  of  the  items  within  the  item  type. 

The  fourth  column  shows  the  number  of  demands  observed  for  each  of  the 
items  in  the  item  type  over  a  six-month  period,  except  in  the  case  where 


ITEM 

TYPE 


NO.  OF 

UNIT 

DEMANDS 

RESPONSE 

ITEM 

NO*  OF 

UNIT 

DEMANDS 

RESPONSE 

ITEMS 

cost  n> 

PER  YEAR 

TIME  (DAYS! 

TYPE 

1  TEHS 

cost  m 

PER  YEAR 

TIME  1  DAYS  1 

1 

2330- 

0*5 

6. 

49 

4 

352* 

7,  f> 

8 , 

3 

1  38  54. 

0*5 

8  * 

50 

1 

900. 

7-0 

8* 

1 

1261, 

1*0 

fl* 

Si 

2 

2  083- 

7-0 

0. 

1 

3  300' 

2*0 

3* 

52 

1 

6200* 

7.0 

8* 

1 

1 7S 1  ■ 

5*0 

8  - 

53 

1 

71* 

8,0 

B. 

4 

01- 

0-5 

0* 

54 

t 

537* 

0,0 

0  * 

7 

214* 

0*5 

8* 

55 

1 

47  38* 

8-0 

6- 

6 

350- 

0*5 

8- 

56 

2 

257, 

9.0 

8* 

4 

628* 

0-5 

a- 

57 

1 

569. 

9.0 

0* 

3 

1613- 

0*5 

3  - 

50 

1 

784* 

9.0 

8* 

6 

171* 

1*0 

8  - 

59 

1 

1500. 

9.0 

0- 

1 

462- 

IlO 

a* 

60 

2 

2176- 

9,0 

8* 

1 

863- 

1*0 

a- 

61 

i 

0864* 

9*0 

8  * 

1 

2223- 

1*0 

a* 

62 

1 

411. 

10.0 

8* 

2 

3219. 

1-0 

a* 

63 

i 

863- 

10,0 

a  ■ 

l 

10678, 

1-0 

8- 

64 

1 

3000- 

10*0 

a* 

1 

196, 

2*0 

8- 

65 

l 

3296* 

10*0 

a. 

2 

388  - 

2-0 

a* 

66 

i 

4800- 

10,0 

8- 

2 

608. 

2-0 

a- 

67 

l 

1404. 

11*0 

8* 

2 

1132* 

2*0 

a- 

68 

1 

1750. 

11*0 

8, 

1 

2  900. 

2*0 

8* 

69 

1 

8500. 

11*0 

8- 

2 

1414* 

2*0 

a* 

70 

i 

379, 

12*0 

8- 

2 

73, 

3-0 

e. 

71 

l 

1297, 

12*0 

0* 

1 

136* 

3-0 

s. 

72 

i 

3007* 

12*0 

8* 

5 

344* 

3*0 

a* 

73 

2 

1240, 

13*0 

8. 

2 

631- 

3.0 

8- 

74 

1 

4400- 

13-0 

8, 

1 

950. 

3.0 

3. 

75 

2 

376. 

14.0 

8- 

2 

1289- 

3-0 

a  * 

76 

1 

1090* 

15.0 

0. 

1 

2147* 

3.0 

8. 

77 

1 

2370* 

15-0 

8  ■ 

2 

2  862, 

3*0 

a  - 

78 

1 

659* 

16*0 

8. 

2 

SO. 

4.0 

a . 

79 

1 

2066* 

16-0 

8  * 

2 

215* 

4,0 

a. 

80 

1 

1246* 

17.0 

8. 

3 

377* 

4.0 

8* 

81 

1 

3898. 

17.0 

B. 

l 

784- 

4.0 

e. 

82 

1 

36* 

10*0 

6  - 

2 

1393, 

4-0 

3* 

83 

L 

796* 

18,0 

8. 

1 

3051- 

4*0 

0- 

84 

1 

112  W 

18*0 

8. 

1 

7458- 

4*0 

a- 

85 

1 

2053* 

10*0 

B  * 

2 

172- 

5-0 

0  * 

86 

1 

675* 

19-0 

a  ■ 

1 

449, 

5-0 

3  * 

67 

1 

3625* 

19-0 

e. 

l 

750, 

5*0 

8* 

88 

1 

9640* 

19-0 

8* 

1 

923* 

5-0 

8* 

89 

1 

13730* 

19,  0 

a  * 

2 

1208* 

5*0 

a* 

90 

1 

132* 

23.0 

a. 

1 

2BB4. 

5.0 

a  * 

91 

1 

700* 

24-0 

8, 

1 

25. 

6*0 

8* 

92 

L 

410* 

27-0 

8  * 

t 

139. 

6*0 

0* 

93 

1 

17600* 

27,0 

8- 

l 

896  * 

6.0 

3. 

94 

1 

47  5* 

31.0 

3- 

1 

1449. 

6*0 

8* 

95 

l 

12570* 

31-0 

a* 

1 

2632* 

6*0 

0. 

96 

1 

124, 

34*0 

8. 

ITEH 

WO.  OF 

UNIT 

DEMANDS 

RESPONSE 

TYPE 

1  TEHS 

COST  [SI 

PER  YEAR 

TIME  (DAYS) 

97 

1 

1830. 

35.0 

e. 

98 

1 

2775. 

46-0 

a. 

99 

1 

1560. 

49-0 

a. 

100 

I 

3430* 

49*0 

a. 

101 

l 

196. 

49-0 

0- 

102 

25 

-  66* 

49.0 

8, 

103 

40 

177. 

49-0 

a. 

104 

30 

355. 

49,  0 

8. 

105 

19 

609. 

49*0 

6. 

106 

2 

863* 

49.0 

8. 

10T 

6 

1174, 

49.  0 

a. 

100 

1 

1568. 

49-0 

a. 

109 

2 

3537. 

49-0 

8. 

110 

13 

62. 

0.5 

30- 

111 

15 

159. 

0.  5 

30- 

112 

15 

342. 

0.5 

30- 

113 

6 

594, 

0-5 

30. 

114 

4 

821. 

0-5 

30* 

115 

3 

1339. 

0.  5 

30. 

116 

1 

2223. 

0.5 

30. 

117 

1 

8143- 

0.5 

30* 

tie 

1 

79. 

1.0 

30- 

119 

15 

157, 

1*0 

30. 

120 

10 

366. 

no 

30, 

121 

6 

559, 

no 

30. 

122 

1 

800, 

1.0 

30, 

123 

1 

1080. 

1.0 

30, 

124 

2 

1703. 

no 

30, 

125 

1 

3056, 

1-0 

30- 

126 

3 

75. 

2 .  Q 

BO. 

127 

7 

215. 

2.0 

30. 

128 

6 

326. 

2-0 

30. 

129 

2 

552. 

2.0 

30, 

130 

1 

855* 

2-0 

30, 

131 

1 

2980. 

2.0 

30. 

132 

2 

40. 

3.0 

30. 

133 

4 

178, 

3-0 

30- 

134 

5 

433* 

3.0 

BO* 

135 

1 

555- 

3.0 

30. 

136 

1 

814, 

3-0 

so- 

137 

1 

1150. 

3,0 

so. 

136 

1 

23* 

4.0 

30, 

139 

2 

163- 

4.0 

30* 

140 

2 

326, 

4-0 

30. 

141 

2 

704- 

4-0 

30. 

142 

3 

1572. 

4-0 

30. 

143 

3 

38- 

4.0 

30. 

144 

2 

217, 

5.0 

30. 

ITEM 

NO*  OF 

UNIT 

DEMANDS 

RESPONSE 

TYPE 

t  TEM5 

COST  IS) 

PER  YEAR 

TIME  ( DAYS  1 

145 

3 

396. 

5.0 

30- 

146 

2 

806, 

5-0 

BO. 

147 

1 

53, 

6.0 

30, 

148 

3 

160. 

6.0 

30, 

149 

i 

345. 

6,0 

30. 

150 

l 

593. 

6.0 

30- 

151 

3 

1844, 

5-0 

30. 

152 

1 

225, 

7.0 

30, 

153 

1 

723. 

7*0 

30* 

154 

1 

1409, 

7.0 

30- 

155 

l 

13. 

8-0 

30, 

156 

1 

283, 

8-0 

30. 

157 

1 

3019. 

8.0 

30- 

158 

2 

72. 

9.0 

30, 

159 

1 

229* 

9.0 

30, 

160 

1 

301. 

9.0 

30. 

161 

1 

849. 

9.0 

30, 

162 

1 

2015- 

9.0 

30. 

163 

1 

3229. 

9,0 

30, 

164 

1 

44, 

10,0 

30. 

165 

2 

175, 

10.0 

30. 

166 

1 

414. 

10.0 

30. 

167 

2 

426. 

mo 

30. 

168 

2 

643* 

mo 

30* 

169 

1 

2266. 

mo 

30, 

170 

1 

32. 

12.0 

30. 

171 

1 

472, 

13*0 

30, 

172 

1 

800. 

13,0 

30- 

173 

1 

371- 

14,0 

BO, 

174 

1 

33, 

15-0 

30, 

175 

1 

1590- 

15.0 

30. 

176 

1 

34, 

16. 0 

30, 

177 

1 

200. 

16,0 

30- 

178 

1 

2769, 

16.0 

30, 

179 

1 

5270, 

16-0 

30. 

180 

1 

512- 

17.0 

30. 

181 

l 

725- 

19,0 

30* 

182 

1 

1457, 

20.  0 

30* 

183 

1 

28500, 

22,0 

BO¬ 

184 

1 

520- 

23-  0 

BO. 

185 

1 

71, 

25.0 

30. 

186 

X 

927, 

26-0 

30. 

187 

1 

1980- 

26.0  v 

30. 

18S 

1 

21B. 

28*0 

30. 

189 

1 

12* 

48.0 

30, 

there  were  no  demands  in  which  case  the  number  .5  appears.  For  each  of 
the  items  within  an  item  type,  the  demand  was  assumed  to  be  Poisson  distri¬ 
buted  with  an  average  rate  per  six  months  equal  to  the  number  shown  in 
this  column.  Each  item  was  assumed  to  have  one  of  two  response  times, 
eight  days  for  items  with  an  ERC  (expendability ,  recoverability,  and 
cost)  code  of  XD1 ,  and  thirty  days  for  all  the  others. 

OBSERVATIONS 

The  results  of  computation  are  summarized  in  Table  2  and  depicted 
in  Figs.  1  and  2.  Inspection  of  the  table  and  figures  leads  to  the  fol¬ 
lowing  general  observations. 

A.  The  NORS  function  is  almost  linear  for  a  wide  range  of 
cost .  This  means  that  a  small  change  in  the  Lagrange 
multiplier  will  result  in  a  markedly  different  stockage 
policy . 

Since  an  economic  interpretation  of  the  Lagrange  multiplier  in  the  con¬ 
text  of  the  NORS  optimization  problem  is  the  incremental  investment  re¬ 
quired  to  reduce  NORS  aircraft  by  one  unit,  one  practical  implication  of 
the  above  observation  is  as  follows: 

The  assumptions  underlying  our  prediction  of  expected  NORS  are  cer¬ 
tainly  not  entirely  true--for  example,  we  cannot  be  sure  that  all  items 
essential  to  the  operation  of  the  aircraft  have  been  included  in  the 
computation.  Our  prediction  of  expected  NORS  is  therefore  probably  not 
accurate.  This  fact,  coupled  with  the  near-1  ineari ty  of  expected  NORS 
as  a  function  of  cost,  means  that  one  should  probably  not  determine  a 
stockage  policy  by  equating  the  incremental  stockage  cost  of  reducing 
expected  NORS  by  one  unit  with  some  estimated  cost  of  a  NORS  aircraft. 


TABLE  2 


CHARACTERISTICS  OF  VARIOUS  STOCKAGE  POLICIES  FOR  FOUR  OPTIMIZATION  CRITERIA 

t INVESTMENT  IN  $  THOUSAND) 


NORS 


INVENTORY 


NVESTMENT 

NORS 

RANGE 

250* 

6*96 

0.40 

500* 

6*39 

0,50 

625- 

6-  19 

0.59 

750* 

6.03 

0-65 

87  5* 

5*83 

0.79 

1000- 

5-73 

0.85 

1250- 

5*99 

0  98 

1500. 

5*22 

0,99 

1750- 

5*03 

0*99 

2000- 

4-93 

1*00 

OPERATIONAL  RATE 


INVENTORY 


INVESTMENT 

NORS 

RANGE 

250. 

7*26 

0.60 

500. 

6.83 

0*86 

615- 

6*76 

0-90 

764* 

6*40 

0*94 

067. 

6*32 

0*97 

999- 

6*02 

0-93 

12  54- 

5*54 

0,99 

1499- 

5*25 

0-99 

1755. 

5*03 

1*00 

1999- 

4*93 

1*00 

BACKORDERS 


INVENTORY 


INVESTMENT 

NORS 

RANGE 

251. 

7-75 

0-57 

348. 

7-73 

0-72 

626. 

6*00 

0-83 

743. 

6*86 

0*80 

074.. 

6*60 

0*92 

1016* 

6*62 

0*95 

1179. 

6-60 

0-97 

1502. 

5.26 

0.98 

1764* 

5-05 

0,99 

2000* 

4*95 

0*99 

FILLS 


I NVENTDftY 


NVESTMENT 

NORS 

RANGE 

249- 

7.83 

0-43 

497  » 

7.48 

0.72 

625. 

7*07 

0.00 

751- 

7, 01 

0,89 

095. 

6,03 

0.93 

999- 

6*32 

0.94 

1249- 

6  .  56 

0.98 

1649- 

5.  15 

0*99 

1750. 

5  *  09 

0,99 

2000* 

4*  94 

1,00 

FILL 

RATE 


0-78 

0-90 

0-93 

0-95 

0*96 

0.97 

0-98 

0.99 

0,99 

0,99 


0  250  500  750  1000  1250  1500  1750  2000 

Inventory  investment  {$1000} 


Fig.  1--N0RS  vs.  cost  of  various  stockage  policies 


Range 


Fig.  2 - - Range  vs.  cose  of  various  stockege  policies 


A  better  procedure  is  to  set  a  constraint  on  the  cost  of  the  policy,  and 
minimize  expected  NORS  aircraft  subject  to  this  constraint. 


B.  In  terms  of  the  expected  number  of  NORS  aircraft,  the  NORS 
policy  is  superior  to  the  other  three,  as  it  should  be. 

The  operational  rate  policy  seems  to  be  the  best  among  the 
remaining  three ,  and  the  fill  rate  policy  is  the  least 
satisfactory . 

C.  For  a  relatively  high  stockage  investment,  the  four  poli¬ 
cies  yield  about  the  same  number  of  expected  NORS  aircraft. 
However,  when  the  constraint  on  investment  becomes  more 
stringent,  the  differences  in  the  characteristics  of  the 
four  policies  ( measured  in  terms  of  the  predicted  NORS) 
become  more  pronounced,  i.e.,  the  performance  of  the 
three  non-NORS  policies  degrades  rapidly  as  the  invest¬ 
ment  in  inventory  decreases.  This  degradation  shows  up 
well  before  investment  has  dropped  below  the  level  needed 
to  support  a  90-percent  fill  rate. 

D.  The  range  of  the  NORS  policy  is  less  than  those  of  the 
other  three  policies  for  a  relatively  low  investment ; 
however,  if  there  is  a  sufficient  investment  fund,  the 
NORS  policy  expands  its  range  more  quickly  than  the 
others . 

The  reason  for  the  above  observation  is  that  when  computing  stock  levels 
based  on  NORS  as  a  criterion,  and  when  available  funds  are  low,  we  de¬ 
pend  on  cannibalization  to  satisfy  demands  for  low- demand  items  and  stock 
more  deeply  the  high-cost,  high-demand  items. 

The  three  non-NORS  policies  discussed  so  far  were  computed  under  the 
assumption  that  no  cannibalization  will  take  place,  We  also  computed 
three  policies  by  optimizing  operational  rate,  backorders,  and  fill  rate, 
respectively,  under  the  assumption  that  one  aircraft  was  available  for 
cannibalization.  Table  3  presents  the  cost  and  performance  characteris¬ 
tics  of  these  policies. 


A  comparison  of  Tables  3  and  1  shows  that  when  investment  is  low* 
each  of  the  three  policies  computed  under  the  assumption  that  one  air¬ 
craft  is  available  for  cannibalization  performs  better  than  the  corre¬ 
sponding  policy  computed  for  the  case  of  no  cannibalizations*  But 
conversely,  when  investment  is  high,  the  no-cannibalization  policies 
outperform  their  corresponding  one-cannibalization  policies.  This 
phenomenon  is  to  be  expected:  when  investment  is  low,  one  should  be 
resigned  to  having  at  least  one  NORS  aircraft  that  can  be  cannibalized. 
When  investment  is  high,  one  should  not  always  count  on  having  a  NORS 
aircraft  for  cannibalization. 


TABLE  3 

CHARACTERISTICS  OF  STOCKAGE  POLICIES  DERIVED  WITH  ONE  CANNIBALIZATION 
( INVESTMENT  IN  $  THOUSAND! 


OPERATIONAL  RATE 


INVENTORY 


INVESTMENT 

NORS 

RANGE 

249. 

7.12 

0.  50 

501- 

6,55 

0.69 

633- 

6,30 

0.77 

7SQ. 

6. 16 

0,79 

876, 

5.93 

0*  84 

999, 

5,75 

0,  89 

1249* 

5.49 

0.93 

1501, 

5,32 

0,97 

1749, 

5,19 

0,99 

1995, 

5.15 

1*00 

BACKORDERS 


INVENTORY 

INVESTMENT 

NORS 

RANGE 

242  - 

7,54 

0,46 

500. 

6.76 

0,64 

62  5, 

6-71 

0,70 

754. 

6.66 

0-78 

876. 

6-61 

0.84 

893  * 

6-60 

0-34 

1250. 

5,50 

0,90 

1505, 

5.34 

0,95 

1750, 

5.21 

0,98 

2002* 

5.16 

0,99 

FILLS 


inventory 


INVESTMENT 

NORS 

RANGE 

253. 

7-72 

0*44 

500. 

7,03 

0.59 

625. 

6-87 

0,68 

751. 

6*65 

0.72 

869, 

6.61 

0*76 

1000. 

6.58 

0-84 

1250. 

5,53 

0.  89 

1475* 

5,38 

0.94 

1749. 

5.24 

0.  98 

2000. 

5*  16 

0.99 

VI .  CONCLUSIONS 


Our  experience  with  the  mathematical  prediction  and  optimization  of 
NORS  is  still  quite  limited,  so  the  only  conclusions  we  can  draw  as  yet 
must  be  tentative.  With  this  qualification  in  mind,  the  following  are 
some  conclusions  that  can  be  drawn  from  the  study. 

1.  Although  we  do  not  have,  as  yet,  an  algorithm  for  optimizing 
NORS  that  la  guaranteed  to  work  in  every  case,  the  algorithm  we  used  in 
our  numerical  examples  seems  to  have  worked  in  every  case  that  we  tried. 
Moreover,  the  algorithm  is  capable  of  starting  with  a  policy  determined 
by  optimizing  some  other  criterion  (backorder)  and  converging  to  another 
policy  which  will  almost  certainly  be  better  as  judged  by  NORS,  and  will 
never  be  worse. 

2.  At  some  levels  of  Investment  (low  in  the  case  we  have  examined, 
but  still  corresponding  to  fill  rates  of  more  than  90  percent)  one  can 
make  appreciable  savings  in  cost  at  no  sacrifice  in  performance  (as  mea¬ 
sured  by  expected  NORS)  by  optimizing  expected  NORS  aircraft  rather  than 
fill  rate  or  backorders. 

3.  Of  the  three  separable  measures  of  performance  (fill  rate,  back¬ 
orders,  operational  rate),  operational  rate  seems  to  be  the  best— as 
Judged  by  expected  NORS. 

4.  One  should  not  attempt  to  determine  a  stockage  policy  by  mini¬ 
mizing  a  total  cost  function  In  which  the  cost  of  a  NORS  aircraft  is  made 
explicit.  Such  a  procedure  is  too  sensitive  to  the  supposed  cost  of  the 
NORS  aircraft.  Instead,  one  should  determine  a  stockage  policy  by  set¬ 
ting  a  constraint  on  investment  and  optimizing  supply  performance,  or  by 
minimizing  cost  subject  to  a  constraint  on  supply  performance. 


Appendix  A 
PROOF  OF  THEOREM  2 


First  let  mb  simplify  notation.  For  each  q  =*  (q^ ,  - . •  ,qn)  €  Q  and 
each  k=0 ,1 ,2 1 . . .  *K,  write 


n 

cpk(q)  »  n  FjCq^+kaj)  . 
J 


Then  theorem  2  may  be  restated  as  follows. 

THEOREM  2*.  Suppose  q*  €  Q  and  let  =  cp^(q*)  for  k=0 . K . 

Then  for  any  other  q  €  Q  we  have 


K  K 

E  £f>k(<0  -  XC(q)  >  £  9k(q*>  -  XC(q*) 
k-0  k-0 


whenever 


K  K 

E  b  log  cp  (q)  -  XC(q)  >  t  log  q>k(q*)  -  XC(q*) . 
k-0  k-0 


Before  proving  the  theorem  we  need  a  lemma. 

LEMMA.  For  any  positive  numbers  t  and  t*  we  have 


t  >  t*[log(t)  -  log(t*>  +  l]. 


Proof  of  Lenma .  Expand  log  t  about  t*  by  a  Taylor  expansion: 


log(t)  -  log(t*>  +  (t-t*)/t*  -  (t-t*)2/2t2, 


-  2-2 
where  t  is  between  t*  and  t.  Since  (t-t*)  /2t  i  0,  the  above  equation 

yields  log(t)  <  log(t*)  +  (t-t*)/t*,  from  which  we  have  our  result. 

Q.E.D. 

Proof  of  theorem.  We  suppose  that  q*  €  Q,  and  bg,...,bK  are  defined 
as  in  the  statement  of  the  theorem*  Also  let  q  €  Q  and  suppose  that 

(i)  Zbfc  log  fck(q)  “  XC(q)  >  Ebk  log  tpk(q*)  -  XC(q*). 

k  k 

We  must  show  that 


E<p  (q)  ~  XC(q)  >  E<pk(q*)  -  XC(q*)* 
k  k 


To  do  this  we  make  use  of  our  lemma  with  cok(q)  in  place  of  t  and  cpk(qA) 
in  place  of  t*  to  write 


(ii)  ^ (q>  >  <Pk(q*)[log  cpfc(q)  "  log  9k(q*)  +  1] 


for  k«0,...,K.  Using  (ii)  first,  then  the  definition  of  b^f.^b^,  and 
finally  (i),  we  have 

EC.  (q)  -  XC(q) 
k  K 

>  Ecpk(q*)[log  «pk(q)  -  log  q>k(q*)  +  1]  -  AC(q) 
k 


=  Ebklog  <Pk(q)  ~  XC(q)  -  Ebklog  <Dk(q*>  +  Ecpk(q*) 
k  k  k 


>  Ebklog  ^(q*)  -  XC(q*)  -  Ebklog  cp^q*)  +  E(pk(q*) 
k  k  k 


fcp^q*)  -  AC(q*) 


Q.E.D 


Appendix  B 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

I.  PROGRAM  DESCRIPTION 

The  program  is  written  in  FORTRAN  IV  G  for  the  RAND  IBM  360/65  in¬ 
stallation,  Because  of  the  limited  accuracy  of  the  IBM  360/65  many  vari¬ 
ables  are  doub le- prec  is  ion . 

Item  data  may  be  aggregated  into  cells  containing  items  with  simi¬ 
lar  cost,  demand,  and  response  time  characteristics.  All  items  in  a 
cell  will  have  identical  stock  levels.  If  optimized  item  stock  levels 
are  to  be  written  on  reserve  tape,  the  user  must  supply  the  appropriate 
Job  Control  Language* 

Current  array  dimensions  allow  up  to  500  cells,  maximum  item  stock 
levels  of  9,  and  up  to  9  aircraft  available  for  cannibalization. 

Approximately  94K  bytes  of  core  storage  are  required  for  execution. 
The  run  described  previously  required  149  seconds  execution  time  for  488 
items  aggregated  into  189  cells,  and  10  targets. 

The  program  consists  of  a  main  routine  and  three  subroutines.  The 
main  routine  handles  input  and  final  output,  and  for  each  target  performs 
four  optimizations  using  as  performance  measures:  fill  rate,  backorders, 
operational  ready  rate,  and  NORS ,  respectively* 

Subroutine  TABLE  computes  the  utility  function  table  for  each  of 
the  four  criteria*  Subroutine  INITL  computes  cost  and  performance  asso¬ 
ciated  with  the  initial  end-points  at  the  beginning  of  each  optimization* 
Subroutine  STEF!  performs  sub-optimization  for  each  cell  during  each  iter 
ation  of  the  optimization  process* 


II.  DATA  DECK  OUTLINE 


The  data  deck  outline  given  below  la  discussed  In  greater  detail 
in  the  next  section. 


FORMAT 

COLUMN 

VARIABLE 

DESCRIPTION 

Card 

1 

15 

1-5 

NTARG 

Number  of  targets 

Card 

2 

15 

1-5 

ITPR 

Print  control 

15 

6-10 

IPSL 

Print  control 

15 

11-15 

ITAPE 

Tape  control 

15 

16-20 

NCELL 

Number  of  cells  of  item  data 

15 

21-25 

KMAX 

Maximum  number  of  cannibal izat ions 

15 

26-30 

MAXSL 

Maximum  item  stock 

F10.0 

31-40 

B 

Length  of  data  period 

Card 

3 

13 

3-5 

NC 

Number  of  items  in  cell 

F10.0 

6-15 

CST 

Unit  cost 

F10.0 

16-25 

DM 

Observed  past  demand 

F10.0 

26-35 

RTM 

Response  time 

Card  3  appears  once 

for  each 

cell  of  aggregated  item  data,  i,e. 

NCELL 

times . 

Card 

4 

10F8.0 

1-80 

A(K)  for 

K-1,10 

Weights  for  the  objective  function 

Card 

5 

F8.0 

1-10 

CONSTR 

Cost  target 

Cards  4  and  5  appear  for  each  target*  i.e.  ,  NTARG  times. 


III.  DESCRIPTION  OF  INPUT  DATA 


Card  I 

NTARG  Any  number  of  cost  targets  can  be  submitted  in  a  single 
run,  hence  NTARG  may  be  any  positive  integer.  The  para¬ 
meters  on  Card  2  and  the  item  data  will  be  the  same  for 
all  targets.  Execution  time  is  proportional  to  NTARG. 

Card  2 

ITPR  If  ITPR  >  0,  output  from  subroutine  TABLE  will  be  printed. 
Output  consists  of  the  demand  distribution,  cumulative 
demand  distribution,  and  utility  function  for  each  cell, 
for  each  of  the  four  optimization  times.  If  ITPR  =  0, 
this  output  is  suppressed. 

IPSL  If  IPSL  >  0,  the  main  routine  will  print  cell  number, 

number  of  items  in  cell,  unit  cost,  past  demand,  response 
time,  and  optimal  stock  level  for  each  cell,  for  all  four 
optimization  types.  If  IPSL  »  0,  this  output  is  suppressed. 

ITAPE  If  ITAPE  >  0,  optimized  stock  levels  are  written  on  data 
set  reference  number  ITAPE.  For  ITAPE  a  6,  output  is 
printed,  otherwise  a  reserve  tape  is  written.  For  each 
target  the  following  will  be  written:  cost  constraint 
in  dollars,  and  for  each  cell,  the  item  type  and  optimal 
stock  levels  corresponding  to  fill  rate,  backorders,  op¬ 
erational  ready  rate,  and  NORS  performance  measures.  If 
ITAPE  -  0,  this  output  is  suppressed. 

NCBLL  If  item  data  are  not  aggregated,  NCELL  is  entered  as  the 
number  of  items.  The  program  is  currently  dimensioned 
for  500  aggregated  cells,  or  500  single  items. 

KMAX  The  program  is  dimensioned  for  KMAX  <  9. 

MAXSL  The  program  is  dimensioned  for  MAXSL  <9, 

B  Number  of  days  over  which  past  demand  has  been  observed. 

Card  3 

NC  If  item  data  are  not  aggregated,  NC  ■  1  for  each  cell. 

CST,  EM,  Apply  to  all  items  in  the  cell. 

RTM 

Card  4 


A(K) *  Utility  function  weights  are  chosen  for  the  number  of 
R*l,10  aircraft  available  for  cannibalization  +  1.  For  one 


available  aircraft,  weights  would  be  set:  A(  1)  =  0, 

A(2)  =  1,  A(3)  through  A(10)  blank.  Alternatively, 
setting  weights  A(l)  =*  .5,  A(2)  =»  .5,  A(3)  through  A(10) 
blank,  would  produce  a  utility  function  which  assumes 
equal  probability  for  0  and  1  available  aircraft. 


Card  5 

CONSTR  The  cost  target  is  entered  in  thousands  of  dollars  and 
is  converted  internally  to  dollars. 


IV.  DEFINITION  OF  VARIABLES 


VARIABLE 

DEFINITION 

A(K) 

Relative  importance  of  the  objective  function  when 
there  are  K-l  aircraft  available  for  cannibalization* 

AA(K) 

Weighted  average  of  two  probabilities,  DEXP(WL(K)) 
and  DEXP(WR(K))  where  DEXP(WL(K>)  is  the  probability 
of  having  fewer  than  K  NORS  with  the  left-hand  stock- 
age  policy,  and  DEXP(WR(K))  is  the  corresponding  prob¬ 
ability  for  the  right-hand  stockage  policy.  When  the 
performance  measure  is  NORS,  A(K)  is  set  to  AA(K)  for 
the  next  optimization  pass.  A(K)  corresponds  to  b(K-l) 
in  Sec.  IV. 

B 

Number  of  days  in  the  data  period. 

C 

Cost  associated  with  an  initial  end  point.  Computed 
in  subroutine  INITL. 

CF(IL) 

Cumulative  Poisson  distribution  with  mean  XLAM . 

CONSTR 

Investment  constraint.  Input  in  thousands  of  dollars; 
converted  in  program  to  dollars. 

COST  (I) 

Unit  cost  for  items  in  cell  I. 

CST 

Unit  cost  for  all  items  in  a  cell. 

DEM(I) 

Observed  past  demand  for  items  in  cell  I. 

DM 

Observed  past  demand  for  all  items  in  a  cell. 

F(IL) 

Poisson  density  with  mean  XLAM,  evaluated  at  the  point 
IL-1,  l.e.,  probability  of  observing  IL-1  demands  dur¬ 
ing  the  response  time. 

FRATE 

Fill  rate  for  a  given  stockage  policy. 

I 

Index  used  to  denote  item. 

IPSL 

Print  control.  IPSL  >  0:  print  item  information  in¬ 
cluding  optimized  stock  levels.  IPSL  =0:  don't  print. 

1ST (I ,IUTL) 

Optimal  stock  level  for  items  in  cell  I,  and  optimiza¬ 
tion  type  IUTL. 

ISWI 

Switch  Indicating  stock  level  change  from  previous  iter¬ 
ation  for  any  cell.  ISWI  =0,  no  change  and  optimiza¬ 
tion  is  complete;  ISWI  =1,  at  least  one  stock  level  has 
changed  and  optimization  may  continue. 

ISW2 

ITAPE 

ITPR 

IUTL 

KMAX 

LEVEL(I) 

MAXSL 

NC 

NCELL 
N't  ARC 
NZS 
P(K) 

PS 

RT(X) 

RTM 

SI 

SA 

SAVE(I  j  J) 
SC 

SLO 


Switch  indicating  increase  in  stock  level  for  a  par¬ 
ticular  cell  during  sub-optimization.  ISW2  =  0 ,  no 
increase;  ISW2  =*  1,  stock  level  increase. 

Tape  control.  ITAPE  >  0:  write  stock  levels  for  all 
optimizations  and  all  targets.  ITAPE  <  0:  don't  write 
tape  , 

Print  control.  ITPR  >  0:  print  utility  function  table. 
ITPR  <  0:  don't  print. 

Optimization  type.  IUTL  =*  1,  fill  rate;  IUTL  =  2, 
negative  backorders;  IUTL  =  3,  operational  ready  rate; 
IUTL  *  4,  NORS. 

Number  of  cannibal izations  +  1. 

Current  stock  level  for  items  In  cell  X. 

Maximum  stock  for  any  item. 

Number  of  items  in  a  cell. 

Number  of  cells  of  aggregated  item  data. 

Number  of  Investment  targets. 

Number  of  items  with  positive  stock  level. 

Probability  of  fewer  than  K  NORS  for  performance  mea¬ 
sure  under  evaluation. 

Percentage  of  items  with  positive  stock  level. 

Response  time  for  items  in  cell  1, 

Response  time  for  items  in  a  cell. 

Old  value  for  expected  NORS  when  performance  measure 
for  optimization  is  NORS. 

Expected  NORS. 

Array  used  to  save  values  for  final  output. 

Current  slope  times  cost  of  a  particular  item.  Used 
in  sub-optimization. 

Slope  computed  during  previous  Iteration, 


SLOPE 

T(K) 

IDEM 

TW 

TW1 

Ti 

U(I,J) 

U0(I) 

UTIL(K) 

V 

WL(K) 

WR(K) 

X 

XL 

XLAM 

XNC(I) 

XR 

Y 
YL 


Slope  =  soda  ted  with  current  end-points. 

Utility  for  the  mixed  policy  when  there  are  K-l  air¬ 
craft  available  for  cannibalization. 

Total  observed  demand. 

Weight  associated  with  left-hand  policy. 

Weight  associated  with  right-hand  policy. 

Weighted  utility  of  mixed  policy. 

Change  in  utility  for  an  item  in  cell  I  when  its  stock 
level  changes  from  J-l  to  J. 

Utility  for  an  item  in  cell  X  when  its 'stock  level  is 

0. 

Utility  of  current  policy  when  there  are  K-l  aircraft 
available  for  cannibalization. 

Utility  associated  with  an  end-point.  Computed  in 
subroutine  INITL. 

Utility  of  left-hand  policy  when  there  are  K-l  air¬ 
craft  available  for  cannibalization. 

Utility  of  right-hand  policy  when  there  are  K-l 
aircraft  available  for  cannibalization. 

Cost  of  current  policy. 

Coat  of  policy  associated  with  left-hand  end-point. 
Average  number  of  demands  over  the  response  time. 
Number  of  items  in  cell  I. 

Cost  of  policy  eeaociated  with  right-hand  end  point. 
Utility  of  current  policy. 

Weighted  utility  of  left-hand  policy. 


ono 


***  S  *1  A  P  I:  M  *♦* 

SUPPLY  OPTIMIZATION  AND  PERFORMANCE  EVALUATION  MODEL 


RIAL* 8  UT ILI 10) ,U(500, 1 0)  , UO (  500 ) , A 1 1 0 ) , T (  13)  , 

l  SLOt  SLOPE  »SC»Y<,YL»Y,TW,TWI,WL,WR,S1,SA,P110) 

?  , X , XL »  XR 

3  .XCRtXCRl.CTARG 

COMMON  U0,U»UML,A,T» 

1  CUSnsuO)#DEM(500)*RT(500) »X »  B, TDEM, 

2  LEVEL (500) ,NCELL , K MAX , MAXSl, XNCt 500 ) 

L 

DIMENSION  WL( 10) ,WR( 10) .1ST! 500,4) , AAI 10 » 

1  ( SAVE ( 10» 1 3 ) 

EUUlVAlfcNCE  I MAXSL*  MAXST) 

C 

c 

READ  5  »  N  T  A  RG 

READ  5,  l TPR, IPSL,1TAPE,NCELL,KMAX,MAXSL,B 
PRINT  4tNCfclL, KMAX, MAXSL, B 

4  FORMAT  ( 1HI,3I5,2F10.0> 

5  FORMAT  (6I5,2F10.0) 

C 

NITM  *  0 
THEM  =  0. 

DD  10  I = 1 iNCELL 
READ  7 «  NC*CST,OM*RTM 
7  FORMAT  I 2X, I3.3F10.0) 

IM  =  MOO I  I ,481 

IF  ( IM.EO. i)  PRINT  6 

6  FORMAT  I IHI/1H-,24X, • I TEM  DATA'/// 

1  LH  ,  HX ,  • ITEM  NO.  OF  UNIT  DEMANDS  RESPONSE'/ 

2  III  ,  rtX » •  TYPE  ITEMS  COST  O)  PFR  YEAR  TIME  (DAYS)'/) 
PRINT  B, I ,NC,CST,DM,RTM 

H  FORMAT  I IH  , 9X , 13 , I  6 , F 1 1 . 0 » F9 . I , F i 1. 0 ) 

NITM  =  NITM  ♦  NC 
XNCU  )  =  NC 
COST  €  I )  *  CST 
OEM (  I  )  =  1>M 
Rim  =  RTM 

10  TDEM  =  TDEM  ♦  DM*FLOAT l NC ) 

12  FORMAT  UH  ,15,4013) 

KX  *  KMAX 

C 

C 

DO  268  ITAR&  =1,NTARG 
READ  (5,15)  A 

14  READ  (5,15, END*270)  CONSTR 
CONSTR  *  CONS  TR*1 000. 

ind  *  n 

OPTIMIZE  FOR  ALL  FOUR  PERFORMANCE  MEASURES 

DO  258  IUTL  =1,4 
KMAX  =  MAXST  ♦  1 
IF  ( IUTL. IT. 4)  KMAX  *  KX 


oo  o  o  o  o  o  non 


C  CONSTRUCT  UTILITY  f-UNCTION  T ABLE 

CALL  TABLE! IUTL, ITPR) 

15  FORMAT  < 1QF8.0) 

WRITE  (6,16)  ( A ( K ) , K* 1 » KMAX ) 

16  FORMAT  l IH1 1 1  OF  1 2 . 6) 

SI  =  0. 

00  17  K - 1 ,  KMA  X 

17  SI  *  SI  ♦  A  t  K  > 

00  IB  K= 1 , KMAX 
13  AIK)  -  A(K)/S1 

SI  =  FLOAT ( KMAX )  -  SI 

INITIALIZE  RIGHT  AND  LEFT  END-POINTS,  COMPUTE  SLOPE 

20  CONTINUE 

00  28  1  =  1, NCE LL 
28  LEVEL! I)  *  0 

CALL  INITLtXR»YR,MAXSL) 

00  30  K= 1 , KMAX 
30  WR(K)  =  UTIL(K) 

CALL  I N I TL ( XL , YL  *  0 ) 

00  50  K*  1 »  KMA  X 
50  WL(K)  *  UTIL(K) 

SLOPE  =  ( YR-YL ) /( XR-XL ) 

PRINT  180,  XL, YL, SLOPE, XR,YR 

Y  »  YL 

BEGIN  ITERATION  FOR  THIS  OPTIMIZATION  PASS 
X  »  0. 

100  ISW1  =  0 

00  120  l*l,NCELL 
SC  *  SLOPE*COST(  I  ) 

CALL  STEP1 ( I , ISW1, SC, LEVEL! I )  ) 

120  CONTINUE 
125  FORMAT  (LH  • AAI 3/ ) 

IF  ( ISW1.EQ.0)  GO  TO  200 

REPLACE  ONE  OF  THE  END-POINTS  WITH  NEW  VALUES  OF  X  AND  Y 

Y  =  0 

00  130  K* 1 , KM  AX 
130  Y  *  Y  ♦  UT I L( K ) *A I Kl 

IF  (X.GT.CONSTRI  GO  TO  L50 
C 

XL  *  X 
YL  *  Y 

00  140  K=1,KMAX 
140  WLIKI  «  UTIL(K) 

GO  TO  170 
C 

150  XR  =  X 
YR  »  Y 

DO  160  K-1,KMAX 
160  WR ( K )  =  UTIL(K) 


ono  on  non  non 


C  COMPUTE  NEW  SLOPE »  REPEAT  OPTIMIZATION  UNTIL  SLOPE  DOESN'T  CHANGE 

170  SLO  =  SLOPE 

SLOPE  *  (YR-YD/I  XR-XL) 

PRINT  180.  XL  t  YL  t  SLOPE  » XR»YK 
180  FORMAT  I2H  I F9. 0, • , • ,F9 .4, • I  SLOPE  *  '.1PEU.4, 

1  ?X» ' ( ' »  0PF9. 0t •t,*F9.4*'),l 

XCR  *  X-CONSTR 
XCRl  =  ,001*CON$TR 

IF  ( DABS I XCR )  .LE.  DABSIXCR1!)  GO  TO  200 
IF  ( OARSI SLO- SLOPE!/ IDARSISLOI+DAttSI SLOPE)! 

1  .GT.  .000001!  GO  TO  100 

200  CONTINUE 

OPTIMIZATION  COMPLETE  —  PRINT  RESULTS 
PRINT  205 

205  FORMAT  ( • OOPT I  Ml Z ATION  COMPLETE'//'  POLICY  CHARACTERISTICS'// 1H  t 
1  T18, 'CURRENT' *T38. 'LOW' ,T5l. 'HIGH* ,T66, • MIXED' // ! 

TW  =  ( XR-CONSTRJ/I XR-XL) 

TW1  *  1 .-TW 
DO  210  K* 1 , KMAX 

210  T(K!  *  TW*WL(K>  ♦  TW1*WRIK! 

NORS  EVALUATION 
SA  -  0. 

DO  215  K* 1 .KMAX 
K 1  =  K-l 

PRINT  225*  Kl,UT(L(K) ,WL(K) * WR( K! * T (K ) * A( K) 

225  FORMAT  I 1H  fI5»5H  CANN.5F15.6) 

AA( K )  =  TW*DtXPf WL(K) »  ♦  T W 1 *DEXP ( WR ( K I ) 

IF  I IUTL.LT .4 )  GO  TO  215 
AIK)  =  AA ( K ! 

SA  =  SA  ♦  1 • - AA IK) 

215  CONTINUE 

COMPUTE  EXPECTED  NORS  FOR  OPTIMIZATION  TYPES  it  2»  3 
IF  I IUTL.E0.4 !  GO  TO  220 
KMAX  *  MAXST+l 
CALL  TABLEOtO) 

CALL  INITLICtVtO) 

218  CONTINUF 

DO  219  K* 1  *  KMAX 
PIK)  =  DEXPIUTIL(K)) 

SA  «  SA  ♦  I  1 .-PI K  I ! 

21*  CONTINUE 
220  CONTINUE 

TV  *  T  W*YL  ♦  TWl*YR 
PRINT  230*  Y,YL,YR,TY 
230  FORMAT  UIH  WTO  UTIL  ,4F15.6) 

IF  (IUTL.EQ.4) 

1  PR  I  NT  235*  (AA(K!*K=l*  KMAX) 

235  FORMAT  I 1H0/10F 12 . 7) 

IF  (IUTL.LT. 4)  PRINT  235.  ( Pf K) *K* 1 1 KMAX ) 

PRINT  236. SA 

236  FORMAT  I'OEXPECTED  NORS  =  ',F12.7//1 

FOR  NORS  PERFORMANCE  MFASURE*  REPEAT  OPTIMIZATION  WITH  NFW  VALUES 
FOR  AIK!  UNTIE  THFRE  IS  NO  CHANGE  IN  EXPECTED  NORS 


non  non  non  on 


IF  (  I UTL • L  T . 4 )  GO  TU  250 
IF  ( S L  .E0.  SA)  GJ  TO  250 
SI  *  SA 
GO  TO  20 
250  CONTINUE 

WRITE  16,252)  IUTL.C0NSTR 
252  FORMAT  ( IH-, I  5 , 5X , ' t • , F  10. 2//  ) 

PRINT  STOCK  LEVELS,  RANGE  OF  ITEMS  STOCKED 
NZ  5*0 

DO  255  1*1, NC ELL 
NC  *  XNCU)  +  .5 
IF  IIPSL.GT.  0) 

1  WRITE  16,256)  I ,NC , COS T (II , DEM ( i) , RT { I ) , LEV  EL (I) 
I  STI I , I UTL )  *  LEVEL! I ) 

255  IF  1  LEVEL { I ) • GT.O )  NZS=NZS+NC 

256  FORMAT  C IM  , 2  I  5 , 2F 15 . 2 , F6 . 0, I  5 > 

PS  -  FLOAT (NZS)/FLUAT(NITM) 

WRITE  (6,257)  PS 

257  FORMAT  ('-RANGE  OF  ITEMS  STOCKED  *  ',F5.31 
SAVE  COST,  UTILITY,  RANGE  FOR  FINAL  SUMMARY 
WCST  =  X 

IF  (IUTL.EQ.4)  WCST*CON$TR 
IF  ( IUTL.GT.il  GO  TO  267 
FRATE  *  TW*YL  ♦  TWl*YR 
IND  *  IND  +  1 
SAVE! I TARG, 13 )  -  FRATE 

267  IND  *  IND  ♦  1 
INDX- 14- IND 

SAVE! ITARG, INDX)  *  PS 
IND  *  IND  +  1 
I NDX*14- I ND 
SAVE( ITARG, INDX)  *  SA 
IND  *  IND  +  1 
INOX* 14- IND 

SAVE! ITARG, INOX)  *  WCST/LOOO. 

258  CONTINUE 

WRITE  TAPE 

IF  ( ITAPE  *EQ .0 )  GO  TO  266 
WRITE  (ITAPE, 259)  CONSTR 

259  FORMAT  (F12.2) 

DO  260  I-l.NCELL 

260  WRITE  (ITAPE, 265)  I ,( I  ST! I , I UTL ), I UTL* 1,4 1 

265  FORMAT  (515) 

266  CONTINUE 

268  CONTINUE 

PRINT  SUMMARY  OF  ALL  TARGETS 

270  WRITE  (6,275) 

275  FORMAT  (lHi) 

PRINT  260 

280  FORMAT  UH-,  57X,  *  TABLE  2'//1 
PRINT  282 


282  FCIRK  A 1  (1H  » 2  1 X,  •CHARACTER!  S  T  I C  S  UF  VARIOUS  STOCKAGE  POLICIES  FOR 
l  FOUR  IjPT  I  MI  Z  AT  I  ON  CRITERIA*) 

PRINT  284 

2ft4  FORMAT  ( 1H  , 46X , » ( I NVE S TME NT  IN  1  T HUUS AN D ) * // / ) 

PRINT  286 

286  FORMAT  tlM  , 1 0 X , • NORS * , 1 5 X , » OPER ATI ON AL  RATE » f 1 3 X , • B ACKOROERS • , 

1  21X»  *F ILLS*//) 

PRINT  288 

2 8H  FORMAT  IlH  , 0 X , 4 I ♦ J NVE NTOR V • # 16X ) , • F ! LL • ) 

PRINT  290 

290  FORMAT  t IH  , 8 X ,4 I • 1 NVE S TMENT • , 2 X , • NORS • , 2 X , • R ANGE • , 2X ) , • R A T E* // ) 
DO  305  l=l,NTARG 

305  PRINT  310  *  I  S  AVF  (  I , J ) , J=1 , l 3) 

310  FORMAT  IlH  , 8 X , 4 1 3X , F 5 . 0, 4X , F4 . 2 , IX , F4 .2 , 2X ) , F4 . 2 ) 

PRINT  275 
CALI  tXIT 
END 

C 

C 

SUBROUTINE  STFP1 1  I  •  I  SHI  ,SC«L) 

C 

C  SUB  OPTIMIZATION  ROUTINE 

C 

REAL**  UT  ILI  lQltUl  500, 10)  ,UOI  500),  At  lOl.TUO)  ,0,SC,X 

C 

COMMON  00 1 UiUT I L  »  A  t  T * 

1  COST! 500), DEMI  500) , RTI 500 ) ,X , B, TDEM, 

2  LFVELI500) ,NC£LL, K MAX , MAXSL , XNC I  500) 

EQUIVALENCE  (  MAXSL, MAXST) 

C 

I SW2  =  0 
C 

20  IF  IL.GE. MAXST)  GO  TO  100 
D  *  -SC 

DU  50  K=l,KMAX 
LK  -  L  +  K 

IF  (LK.GT. MAXSL)  GO  TO  50 

D-D+UI l«LK)*A(K) 

50  CONTINUE 

IT  ID.LF.O.)GO  TO  100 
C 

L  =  L  +  l 

53  FORMAT  IlH  ,2110) 

ISWl  =  1 
00  60  K* 1 , KMA  X 
LK  *  L*K-1 

IF  (LK  .GT.  MAXST)  GO  TC  60 

UTILIK)  »UTILIK)  ♦  Ull,LK)  *XNCU) 

60  CONTINUE 
C 

X  *  X  ♦  COSTI 1 )*XNCI I ) 

ISW2  *  l 

GO  TO  20 


100  IF  I  1 SW2  .EU.  1)  GO  TO  200 
120  IF  IL.LF.O)  GO  TO  200 
D  =  SC 

Of)  15(i  K  =  1  ,KMAX 


C 

C 


LK  =  L+K-l 

IF  (LK.GT.MAXSL)  GO  TO  150 

D*D-U< I,LK)*A(KI 

150  CONTINUE 
C 

IF  ( D.LE .0. 1  GO  TO  200 

L  =  L-l 

ISW1  »  1 

00  160  KaUKMAX 

LK  a  L+K 

IF  ILK  .GT.  MAXST)  GO  TO  160 

UTILIK)  *UTIHK)  -  UlItLK) 

160  CONTINUE 
C 

x  =  x  -  cosrm*xNcm 

GO  TO  120 
200  CONTINUE 
C 

RETURN 

END 

C 

C 

SUBROUTINE  I N I TL! C fV, LMM) 

C  INITIALIZE  COST  AND  UTILITY  FUNCTION 
C 

REAL*®  UTILI 10) .UI500, 10) *U0! 500), A! 10),TI 10) » 
1  f C»CTLM, X 

C 

COMMON  UOtUtUT I L  t  A ,T« 

1  COSTI500) , DEMI  500) , ATI  500) ,  X.  B.TOEM, 

2  LEVEL ( 500 1  ,NCELL » KMAX, MAXSL f XNCI 500 ) 
EQUIVALENCE  ( MAXSL* MAXST) 

C 

DO  10  K*1 *KMAX 
10  UTILIK)  «  0 

00  30  I* l .NCELL 

UTIHU  *  UTILI  1 )  ♦  UOU)*XNCm 
LM  «  MAXOI LMMf LEVEL! m 
IF  (LM.EQ.OI  GO  TO  30 
00  20  J»l,LM 

20  UTILtll  *  UTILIll  ♦  Ul t , J l *XNC 1 1 1 
30  CONTINUE 

IF  IKMAX.LE.l)  GO  TO  100 
C 

00  60  Ka2*KMAX 

00  50  I«1,NCELL 

LM  «  MAXOI LMM, LEVEL! I) ) 

KLMl  «  tULM-l 

IF  IKLM1.GT.MAXST)  GO  TO  50 
UTILIK}  *  UTIL(K)  ♦  U( I ,KLM1 ) *XNC( I ) 

50  CONTINUE 

60  UTILIK)  «  UTILIK)  ♦  UTILIK-1) 

C 

100  V  -  0 

DO  120  K*l,KMAX 
120  V  -  V  ♦  A I K I *UT I L I K) 

C  *  0 

DO  150  fal , NCELL 

LM  =  MAXOILMM, LEVEL! I ) I 


♦XNCI I ) 


on  o  on  o  noon 


XLM  =  LM 

CTLM  =  COST! I )*XLM*XNC!I) 

150  C  =  C  ♦  CTLM 
160  CONTINUF 

PRINT  700,  C, V, (UTIl(K) ,K*1,KMAX) 

200  FORMAT  I IH  ,  F  1 0. 0 ,  1  OF  1  2 .5  > 

RETURN 

END 

C 

C 

SUBROUTINE  TABLE! IUTL, I TPR) 

COMPUTE  UTILITY  FUNCTION  TABLE 

REAL*8  UTIL! 10) ,U( 500, 10) ,U0l  500),  A! 10 ) , T ! 10 ) , XLAM, F, CF 
1  ,  X 

COMMON  UOtUtUT IL , A ,T , 

1  COST! 500) , DEM! 500) ,RT! 500) ,X, B,TDEM, 

2  LEVEL  I  500) ,NCELL, KMAX, MA XSL , XNC I  500 ) 

EOUI VALENCE  ! MAX SL , MAX S T) 

DIMENSION  F(20) , CF ( 20 ) 

MAXI  =  MAXST  ♦  1 
Ml  *  MAXST-1 
DO  510  1  * l ,NCELL 
XLAM  =•  DEM! I ) *RT ( I )/8 
F(  1)  =  DEXP! -XLAM ) 

CF! 1)  =  F! 1) 

IF  IMAXl  .LE.  1)  GO  TD  50 
DO  40  I L  =  2 , MA  X l 

FIIL)  *  F( IL-1 )*XLAM/FLOAT ( I L-l ) 

40  CF!  IL)  =  CF!  I  L-l )  ♦  F!  ID 
50  CONTINUE 

IF  ( I TPR  .GT.  0) 

1 PR  I  NT  60,  XL AM , I F ( IL), 1 L-l, MAXI) , (CF! IL) , IL=1,MAX1 ) 

60  FORMAT  < 1 HO, 1 OF  1 2 . 7/ 1H  ,12X,9F12.7) 

GO  TO  ( 100,200,300,300) ,IUTL 

FILL  RATE  OBJECTIVE  FUNCTION 
100  U0( I )  =  0 

DO  120  J*l, MAXST 
120  U(I,J)  =  F(  J)*DEM(  D/TDEM 
GO  TO  400 

200  UO! I )  =  0 

DO  220  J«l, MAXST 
FJ  *  F! J+1)*FL0AT( J) 

U(  I » J )  =  FJ 
22 0  U0(! )  =  UO! I )  -  FJ 
GO  TO  500 

□P.  READY  RATE  OBJECTIVE  FUNCTION 
300  UO! I )  =  DLOGICF! 1 » ) 

DO  320  J*l, MAXST 
320  U(I,J)  =  DLOGICF! J+l) ) 

DO  330  J= 1 , Ml 
JJ  =  MAXST-J+1 

3  3  j  UU,JJ)  =  U  (  I  ,  J  J  )  -  UltfJJ-1) 


ui  i ,  n  =  ui i *  n  -  uo t  r  i 

400  IF  (HAXST  .LF.  1 )  GO  TO  500 
DO  470  J=l,Ml 

IF  (U<  I,J  +  U.LE.UU,  JU  GO  TO  500 
XJ  =  J 

UJ  =  (XJ*U(I,J)  +  U(  I «  J  +  l ) ) / ( X J+l ) 

m  r,  j+n  =  uj 

DO  450  JJM.J 
450  U(I, JJ)  -  UJ 
470  CONTINUE 
500  CONTINUE 

IF  ( ITPR.EU.O)  GO  TO  510 
PRINT  505 i  UQ< I), ( U ( I ,J),J=l,MAXST) 
505  FORMAT  UH  ,  7F  1  6 . 7  ) 

510  CONTINUE 
RETURN 
END 
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